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PREFACE 


This investigation was conduoted at the United States Naval Post~ 
graduate School, Monterey, California as the thesis requirement for the 
degree of Master of Science in Aerology. 

For help and advice received in its preparation the author is indebted 


to Assooiate Professor Ge J. Haltiner of the Postgraduate School. 
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I. INTRODUCTION 


In the attempt to subject the state of the atmosphere to a mathe- 
matical analysis, and in particular, to explain the growth and behavior 
of wave like disturbances in the westerlies, the meteorologist is cone 
fronted with two major problems: 

(a) An imperfeot knowledge of the instantaneous state of the 

atmosphere and 

(vo) The mathematical difficulties inherent in the attempt to 

integrate the equations of motion. 

Adequate observational data is available to permit us to draw some 
conclusions as to the "average" conditions in the atmosphere, or, in the 
case of choosing a particular model to analyze, a "plausible" condition. 

The mathematical problem is simplified by considering the effect of 
a small perturbation on a known zonal flow and neglecting all differentials 
in the resulting equations of higher order than the first. This results in 
a set of linear differential equations <= the perturbation equations. The 
mathematical details have, in most cases, been further reduced by making 
certain simplifying assumptions in the models. Although such assimptions 
introduce artificialities when compared with the atmosphere, they essent- 
Jally resolve the problem to the study of dynamic instability -~ that 
portion of the stability or instability contributed by the nature of the 
flow pattern itself. Kuo [4] states that, in studying two-dimensional 


nondivergent flow, 


(1) 
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Similar results would be obtained for the very general case 

with vertioal motion, solenoids, and divergenoe and oonvergence, 

if certain approximations sre made in reducing the partial dif- 

ferential equations for a certain dependent variable to an ordi- 

nary differential equation. 

The technique of using the perturbation equations is not new, having 
been used over fifty years ago by Helmholtz, Rayleigh and others. In 1928 
Solberg showed that unstable waves oan develop in a sloping surface of dis- 
continuity between two parallel streams of different density, supporting 
the results of his synoptic studies of wave development on the polar fronts 
In 1939 Rossby developed his widely known equation for the speed of long 
waves in the westerlies in the special case of constant zonal motion in a 
homogeneous incompressible atmosphere. Since that time, a great number of 
meteorologists have applied the theory to a variety of modelse Machta [s] ; 
using Rossby's model, assumed a solution for the perturbation stream func= 
tion which allowed for a tilt to the trough line. Haurwitz also used Rossby's 
model but allowed for the curvature of the earth and obtained only slight 
modification to Rossby's results. Charney (iJ considered a special case of 
baroclinic three-dimensional flowo 

Kuo [4] developed general stability criteria for a two-dimensional, 
barotropic nondivergent atmosphere using general, symmetric velocity profiles. 
Haurwitz |3] considered pure shearing waves between two constant zonal 
two-dimensional wind fields on 4 rotating planet. The oase of a shear zone 
of finite width was considered neglecting the variation of the ooriolis para- 
meter with latitude. This problem in considered briefly in this paper. 

It may be noted that in nearly all of the papers studied only a portion 


of the complete problem was investigated. Some are concerned with the question 


of whether or not a perturbation will be unstable. Others have eliminated 


(2) 
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potentially unstable waves from their solutions (as will be done in this 
paper) and proceed to establish the speed with which the long stable waves 
will move. This approach is entirely reasonable because if the wave is un- 
stable, the solution is valid only in the initial development and the sig- 
nificance of the real values of wave velocity is limited to this time. Also, 
as Charney [1] points out, the evaluation of the complex roots may present 
great difficulties. On the other hand, if the wave is essentially stable, 
the seoond order differentials which have been neglected in the pertur- 
bation equations will continue to be small and the value of wave velocity 
obtains is valid for an appreoiable length of time. 

In this paper, only the long, slow moving or retrograding waves are 
considered. The frequency equation resulting from the solution of the 
boundary systems determined by a sinusoidal velocity profile and a linear 
shear zone are determined. The frequency equation for a probability ourve 
profile was not fully developed due to mathematioal difficulties encowmtered 
in the application of one of the boundary conditions. This matter is dis- 
cussed in some detail. 

The roots of the frequency equation in the case of the sinusoidal 
velocity profile are evaluated (for a particular wave length) for a synoptic 
situation in which the horizontal wind velocity profile was in fair agree- 
ment with the model. The results are found to be of the right order of 
magnitude and show that there may be more than one possible value of wave 


speed under a given set of conditions. 


(3) 
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We shall now consider three velocity profiles and derive the equations 
relating wave speed and wave number. The perturbation equations to be used, 
as developed by Haurwitz [2] , assume incompressible flow. We shall further 
assume two-dimensional nondivergent flow. The variation of the coriolis 
parameter with latitude is assumed to be a linear fumotion of y. (Appen- 
dix II). Under these conditions, the perturbation equations reduce to 


ees) + Vv (0 —=8 Eee ee USO 


(1) 


as derived in Appendix Ie 

A —s analysis of this equation must include a study of the be- 
havior of v in the neighborhood of the singular point c =U. Such a point 
must ocour if the wave is moving with a speed between U, and Ungyx- Kuo (4] 
has investigated this point for the case of general symmetric velocity profiles 
by including friotional terms in the above equation. These terms are small com- 
pared with the others except when o - U is small. 

He shows that if a singular point exists, a stable wave oan exist only 
if U-B=0 at some point on the velocity profile. If this condition is 
met, the speed of such a stable wave is the value of U at that point. Kuo 
calls this the oritical velocity. Any wave moving with less than this 
velocity (but greater than U,) will be unstable, while those moving with a 
velocity greater than the critioal velocity will be damped. He further shows 
that for the long slow moving waves (c 4U,), the amplifying or dampening 


factors are small or zero and such waves are essentially neutral. 


(4) 
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(9) 





In the solutions to be disoussed, no closed form was found and the 
power series solutions obtained are valid only for (o = U) < 0, the slow 


moving or retrograding waves.* 


« Xuo [4] has shown that the minimum wave length of such waves is about 
o200 km under typioal shear oonditions. 


(5) 
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II. SINUSOIDAL VELOCITY PROFILE 


The first problem to be considered is the solution of equation (1) 


where 


+H 


ec eee ° 
(2) 


AQ 
ee ’ 
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+H ° 


A series solution for v in the sinusoidal region is obtained in the oon- 
ventional manner as shown in Appendix III. It involves two arbitrary 
constants which must be eliminated by the imposition of two boundary 
conditions. 

Kuo [4] assumes that for symmetrio velocity profiles and symmetric 
boundary conditions, there will be either symmetry or antisymmetry in the 
amplitude of the perturbation stream function, and therefore only half the 
zone need be considered. He then points out that it would be wmlikely to 
have v = 0 along the line of maximmm wind, thus reducing the study to one 
of symmetric disturbances. 

Kuo's development is in terms of wy , the perturbation stream function 
instead of v. It oan be shown that, wnder the assumptions involved, equa- 
tion (1) in either variable is identical. 

The above statements result in the conclusion that v has a maximum 


where U is a maximum; that is, 
ne 
v (am) = 


(In this and subsequent work, the quantity in parentheses refers to the 


point at which the funstion is to be evaluated) 


(6) 
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This equality is used as the first boundary condition and results in 
the elimination of one of the arbitrary constants. 
Where y < 0, in the field of constant zonal wind, it can be shown 


that the solution of (1) is 


Then 


This condition must be true at all points in the zone, therefore it is 
true at y = 0. Since y = 0 is common to both this zone and the sinusoidal 
zone, the statement is true at this point when applied to the sinusoidal 


zonee The second boundary condition therefore is 


ano) = AV (0) <2 a0 


The imposition of the boundary conditions is shown in Appendix III 
and results in a fourth degree equation ind, where d =¢ - UU. The four 
roots of this equation may be evaluated, by rather tedious methods, after 
selection of the physical parameters involved. 

Two synoptic situations in which the velocity profiles at the SOO mb. 
level were in fairly good agresment with that used in this development are 
considered as examples. These are illustrated on the following page. The 
parameters of the upper profile were evaluated and the equation determined 
for a wave length of 5000 km. This resulted in 

4° - 454 - 1810 4 ~ 240d - 5300 = 0 


The roots of this equation are (approximately) 


ad w <2 m/seoe, ~1) .8 m/S6C., =22 m/seC.e, Tio m/sec. 


(7) 
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The last root may be discarded since we are not considering waves whioh 
move faster than minimim wind velocity. This gives the possible values 


of wave speed as 13.4 m/sec., 13.6 m/sec., and -5.6 m/sec. The wave is 


stable since al] the roots are real. 


(9) 
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III. PROBABILITY CURVE VELOCITY PROFILE 


The second velocity field to be considered in the solution of (1) 


is 
(3) 


By varying the parameters, this profile may be adjusted to give a good fit 
to a large number of symmetric profiles in the vicinity of the jet. The 
differential equation is solved by the series method in Appvendix IV. The 
first boundary condition is the same as in the previous problem and when 
applied results in a solution for v involving one arbitrary constant. 


The second boundary condition, as before, is 
vi = rv = 0. 


However, the "point" at which this must be applied is at infinity. The 
use of this condition therefore requires a transformation of the equation 
in which y == « The point at infinity is then transformed to the origin 
and the boundary condition may be applied. However, the mathematics here 
is rather complicated and the actual expression for c in terms of d was 


not obtained. 


(10) 
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IV. LINEAR SHEAR ZONE VELOCITY PROFILE 


The final problem is conoerned with the solution of (1) for a wind 


profile which is defined by the equations 


veo, 4 Sa 4 of 4 6 L 

(4) 
v= OU, 4 Ss 0 
v= VU. ue 34. 


This model differs from the previous two and from those considered by Xuo 
in that it is not symmetrical. It is also obviously less realistic. The 
solution for v in the shear zone is given in Appendix V. 

By similar reasoning as in the preceeding problems, the bowdary con- 


ditions are found to be 


v @) -rA,vlo) = og 


yo) + navel) = o 


where r) refers to the zone of velooity U); and rj, to the zone of Uj}-> 

By applying these boundary conditions as indicated in Appendix V, a 
rather complicated equation containing an implicit relationship between o 
and & results. The difficulties of finding the roots of this type equation 
are so great that no attempt was made to obtain numerical results for a par- 


ticular model. 


(12) 
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Ve. CONCLUSIONS 


As with any second order differential equation with variable coef- 
ficients, the equation solved in this paper does not, in general, yield 
a solution in terms of elementary fimctions. It is quits possible that 
the problems considered here do possess such solutions but none was founde 
An attempt was mde in each case to transform the equation so as to obtain 
a solution in terms of Bessel or legendre fictions, with no Success. An 
arbitrary selection of wind profile will, in general, be met by the same 
difficulties. 

The perturbation equations are, therefore, readily integrable only 
if the velocity profile is chosen specifically to achieve this purpose. 

A particular but arbitrary profile such as one of those chosen for this 
paper will, in general, yield an equation in which the algebraic diffi- 
culties in obtaining numerical results are of such magnitude as to make 
them practically worthless. These diffioulties would be inoreased in the 
case of complex roots. 

The restriotive assumptions made in all published studies in this 
field have been justified, to a certain extent, but must not be forgotten 
when attempting to apply results to the atmosphere. The processess of the 
real atmosphere are much more complex than we have assumed. The three 
dimensional nature of the atmosphere, divergence, vertical motion, and 
non-homogeneity certainly play an important rols, and often a dominant 


one in the development and propagation of disturbanoes in the zonal flow. 
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APPENDIX I. 


DERIVATION OF SECOND ORDER ORDINARY DIFFERENTIAL EQUATION 
IN v AND y FROM FERTURBATION EQUATIONS 


The perturbation equations for two-dimensional nondivergent flow 


in an incompressible fluid where the undisturbed flow is zonal are: 


du du eU ' o@ 

—— U — sagliaie — $ = _! Sey. 

= * ia 5 v ek (5) 

ov dv - i Sb 

oe oe rte Foy ) 

Mie _ 5 (7) 
x 34 


In order to study the velocity of wave perturbations in the x direction, 


we may assime the following form of solution: 


ca (y-C 2) 
. - “iq 2 (8) 
ce (y-Cct) 
>~v(y4) 2 (8) 
ca (¥-eT) 
(10) 


py) 2 


Then from (7) and (8), 


« dv 
v7) a ory (1) 
Substituting (11) and (10) in (5) we obtain 
fucy-et) 
e(U >) v| (12) 


= eo pene vy) 2 
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Taking the partial derivative of (12) with respect to y and equating it 
to the left hand member of (6), using (9) and (11), we get an equality in 
which each term contains the coefficient ~_" - Dividing through 
by this quantity reduces the equation to one in which the only dependent 
variable is the amplitude factor of v, v(y). For sake of brevity, v(y) 
will now be designated as simply ve The differential equation then 


reduces to 


vi (¢-U) + v\v"- B -a*Ce-v) | ane (1) 


This may be considered as an ordinary differential equation, but it mst 


be remembered that the final solution of v contains the periodicity factor 
cay = CL) 


(15) 
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APPENDIX II. 


VARIATION OF CORIOLIS PARANETER WITH LATITUDE 


Let the coriolis parameter be represented by 


§ = ae + © > 
where $ is assumed constant in the zone wder consideration. Then §G= x 
but 5 = V4 ~~ Dam ® 
(A) 
and as = Gee, Cos ® 0_® 
a4 04 


From the sketch we see that 


| - Kido dy > 
d 
a eo. i 
39 iy 
Substituting this value of = in (A) we have 
cha = B = Ae Zz w co @ 
J R 


(16) 
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APPENDIX III. 
DEVELOPMENT OF FREQUENCY EQUATION FOR SINUSOIDAL VELOCITY PROFILE 


We shall now solve the equation obtained in innate & for the case 


of a sinusoidal velooity profile: 


U= VD, +0, tim by SS hy 
(13) 
The sinusoidal profile is bownded on both sides by a field of constant 
velocity, Upe 
Substituting (13) in (1), we have the following equation to be 
solved: 
dee day) avi Og pam y eh) = G (14) 
where 
dies Cia Ome 
6-4 - - 
es Dien ds 
Expressing sin by in series form 
3 2 
San by. Ly _ (Ly) t (hy) 
34 ay (15) 
and assuming a power series solution for v, 
Vz @o+ Gay ee / mage Ee ee (18) 


We may substitute (15) and (15) in (14) and obtain an equation in positive 
powers of y equal to zero. Upon equating the coefficients of each power 


of y to zero, we find that 


; or 
@,-=-- @, ad 
(-pby-obdy) na 
ne ats lod : Eq 
3 
oe MS) oan), (17) 
THOR V2 4,+4,y-4.44 ¥ - ed - 


(17) 
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This expression for v involves two arbitrary constants which may be elim- 
inated by the imposition of two boundary conditions. The first condition 
to be applied, as shown on page 6, is that v have a maximum where U is a 


@ maximum, or, 
es 
v (aE mae 


Taking the first derivative of (17) with respect to y, setting y ae 


and equating the result to zero, we have 


= Ar Ghee) — Z 
i i eo) ee = 
pig Oe bd > dt yf? (18) 
and upon solving for a, we see that 
um@htd + VU, Rw (h+od) 
6, 2.4, Sek ee =:~SC—t*=<“‘C<;7«C2SC ~*! (19) 


Std -hdT” 
The resulting solution for v, neglecting, for the moment, the third and 


higher powers of y, is 
ythbd +, L(h+oadT ee 
a Ce... ean ve al 


/ = OG Loe oes Laz * a (20) 
By imposing the second boundary condition (see page 7), 
Vato) re {5 } = 6 
we obtain 
a ‘Pu 
uT@hted +ub (heod)e B 
aa  tt*”*Ct—‘i—CC ieee ane = 0 
Ao @2~ 4>— hh dart : d (21) 


This step shows that only the first two terms of the series solution for v 


are Significant in the development of the frequency equation so the powers 


(18 ) 
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of y higher than the first could have been ignored after having deter- 
mined equation (19). However, it should be noted that the solution of 
a, in equation (19) is approximate since only the first four terms of 
the general solution for v were used in its determination. It follows 
that the frequency equation is approximate. 

Squaring equation (21) and collecting the coefficients of powers 
of d we have 

u 
p'(-e Ps ey") (ant cy ; rary 
(22) 
‘i d*( bar'h"s a rig te f a ae grag -1' Pa") 


= a(&rt Bg +27 hq 8 vB) ag 6 = 9 


(19) 
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APPENDIX IV. 


DEVELOPMENT OF FREQUENCY EQUATION FOR PROBABILITY 
CURVE VELOCITY PROFILE 


In this section, the differential equation 


Vi tceOyenNOr =a ene ea) =" ¢ (1) 


is to be solved for the case in which the undisturbed wind velocity 


profile is defined by the probability curve 





ey 
jy = (OR Ui (23) 
A power series solution for v is assumed, 
vz Qo + &, Y tai.y + i a (24) 
and 2. ’ 48 expressed in series forn, 
ee -ky” + Y ae ur a 


The series expressions for v and U are then substituted in (1). The 


4 


coefficients of each power of y are equated to zero and we may now obtain 


expressions for any a, in terms of a, and aj: 
1b) 7 Se e 
eae «a Un, -2RU, ae 6) 
2 (0, ~¢) 
S BUR een “eo i) 
z — where U..-U.+U, and 
aml Cone) 
a, and a, are arbitrary constants. The general solution for v is therefore 


dV, ~ U be = < Ps 
ee aoe La PV 
de oe ee eae 


V=>@,rG@,y~y+e 
vA Q MO} é b (0 
m~ Cc) 


Tae (26) 


(20) 
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We may eliminate one of the constants by the imposition of the boundary 
condition 

Veo) = G6 
The use of this condition for symmetrical velocity profiles is discussed 
on page 6. 
Taking the first derivative of (26) with respect to y, and setting y = 0, 


we obtain 


QA, = 0 


Equation (26) then reduces to 
ene OU Eca auc y” 


] Ea ae 
Se ay] + 2 ( Um-c) 


As mentioned on page 10, the application of the second boundary 


condition, 


v (==?) - nv(-~) = a 
leads to mathematical difficulties and the elimination of a, and consequent 


development of the frequency equation was not accomplished. 


(21) 
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APPENDIX V. 
DEVELOPMENT OF FREQUENCY EQUATION FOR A LINEAR SHEAR ZONE 


In this section the differential equation 


coe ye 2B - Roc atu) 2% (2) 


is solwed for the case in which the undisturbed wind velocity profile is 

a linear shear zone. The arbitrary constants are eliminated and the fre- 
quency equation developed for the case in which the shear tone is termi- 

nated, on each side by a field of constant zonal wind. 


The shear zone is defined by 
Ou -U, 


4 (27) 


b= O, + ik y where Ix = 


A power series solution for v is assumed 
oko Gay + Oy y igi owe 
and this series is substituted in (1). The coefficients of each power 


of y are equated to zero and the following relationships obtained: 


Paap 
Om :- a, 2 


L 
d 








a,-= @ 
ss ° od ~ <F bd 
The series solution for v is then 
Pe PoE, y~ i 8 p d (~*d+3) 3 
= @ Saas +@, See 
Vz @,+ Gy + @, aS: os ae. a. yr (28) 


As explained on page 11, the boundary conditions to be imposed are 


Vi ee) 


os ( £) t 2, y(£) 


i 
>) 


(29) 


iL 
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(30) 
(22) 
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Py substituting (28) in (29) and letting y = O we find that 
a, = we, 2s 


Equation (28) then reduces to 


| | Ko +rid(d'd+B8) a 
\ ei ue 





To impose the second boundary condition, we substitute (31) in (30) and 


set yee After dividing through by the common factor ay, we obtain the 


frequency equation which becomes, upon rearranging: 


bd (R+2,4nn,$) + ba} (ids B) 


t (32 
i ine) ae lee sD + ad'f ween (4 A+ a). | 


t Ny Deer, d(atde 4) | =o 
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This is approximate since fourth and higher powers of y have been neglected. 


(23) 
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